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Abstract. We report on the search for the Standard Model Higgs boson decaying into two photons based on 
proton-proton colhsion data with a center-of-mass energy of 7 TeV recorded by the ATLAS experiment at the 
LHC. The dataset has an integrated luminosity of about 1.08 fb^'- The expected cross section exclusion at 95% 
confidence level varies between 2.0 and 5.8 times the Standard Model cross section over the diphoton mass range 
110 - 150 GeV. The maximum deviations from the background-only expectation are consistent with statistical 
fluctuations. 



1 Introduction 

In the Standard Model (SM), the introduction of a scalar 
field with nonzero average value is responsible for the spon- 
taneous Electro Weak Symmetry Breaking (EWSB), giving 
masses to the intermediate W"^ and bosons and regular- 
izing the weak bosons quartic coupling divergence at high 
energies This scalar field leads to the theoretical ex- 
istence of a further elementary massive quantum particle 
called the Higgs boson. In spite of efforts of past experi- 
ments at LEP and the TeVatron, the Higgs boson has not 
yet been observed |2|,|3|. Several mass ranges have how- 
ever been excluded at 95% Confidence Level, including 
Higgs mass regions below 1 14.4 GeV and from 141 to 476 
GeV 0,11. 



2 SM Higgs boson decaying to two 
photons 

Following the Standard Model predictions, the Higgs bo- 
son can be produced by the LHC proton-proton collision 
by several processes. The main ones are: the gluon fusion 
through a quark loop, the vector boson fusion, the asso- 
ciate production with a or Z" and the associate pro- 
duction with a top-antitop pair. The corresponding produc- 
tion cross-sections as a function of mass are documented 
in Ref. |i6J. The SM Higgs boson can decay to a multitude 
of final states, including a pair of bottom quarks, tau lep- 
tons or intermediate bosons, etc. Despite its small branch- 
ing fraction of 0(1Q^^), the diphoton final state has the ad- 
vantage of allowing the full reconstruction of the Higgs 
candidates mass while keeping the background to an ac- 
ceptable level. 



3 IHiggs candidates selection 

The results presented here can be found in Ref. |7|. The 
ATLAS detector is used to reconstruct pp collision events 
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given by the LHC. They are selected by a primary dipho- 
ton trigger requiring two photon objects with a transverse 
energy Ej greater than 20 GeV. Events are required to have 
at least two offline reconstructed photons. The final dataset 
used in the diphoton analysis represents a total integrated 
luminosity of approximately 1.08 fb^^- 

The event selection relies on several ATLAS subdetec- 
tors: the pixel detector, the semiconductor tracker, the tran- 
sition radiation tracker and the electromagnetic calorime- 
ter where the photons deposit most of their energies. Both 
converted and unconverted photons may enter the dataset. 
One event must have at least one primary vertex with 3 
associated tracks, for which the individual transverse mo- 
mentum pj is greater than 0.4 GeV. The photons, recon- 
structed from the electromagnetic clusters in the event, are 
required to be outside of a region I77I € [1.37, 1.52] which 
represents the transition region between the electromag- 
netic calorimeter's barrel and endcaps, and inside |//| < 
2.37. The photon with the highest transverse momentum 
is called the leading photon, while the one with the sec- 
ond highest is called the subleading photon. The first is re- 
quired to have a transverse momentum p^'^'^ ~^ greater than 

40 GeV while the second pij,"''''^'*'' has to be greater than 
25 GeV. The photons are asked to be isolated. The iso- 
lation is computed from the energy contained in a (77, (p) 
circle around the photon candidate, the core region of the 
electromagnetic cluster being excluded from the calcula- 
tion. For each photon, it is required that El^° < 5 GeV. 
Both photons are asked to pass tight selection criteria, 
which relie on shower shape variables and energy leakage 
inside the hadronic calorimeter. These identification crite- 
ria have an efficiency of about 80% for true photons. Fi- 
nally, only the events with 100 < niyy < 160 GeV enter 
the final dataset. Events passing this offline selection are 
selected by the trigger with an efficiency greater than 99%. 



The final ATLAS dataset consists in 5063 events, for 
which the diphoton invariant mass niyy is plotted on Fig. 
[1] The background is compared to the dominant yy, y- 
jet and jet-jet processes MC predictions. The background 
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level and rejection compare well to the QCD predictions 
and detector performances fSl. To study the SM Higgs sig- 
nal, several Monte-Carlo (MC) generators are used. The 
PowHeg is used to simulate the signal events from the 
dominant gluon fusion process. MC @NLO is used for cross- 
checks. PowHeg is also used for the vector boson fusion 
process, while PYTHIA is used for associate production 
processes with a vector boson or a top-antitop pair A full 
simulation of the detector geometry and response with the 
GEANT 4 program is performed. The samples are pro- 
duced so that the multiple interaction pileup is close to 
the one observed in proton-proton collision data. The same 
selection is applied on MC and data samples. The gener- 
ator cross-sections are normalized to the NNLO predic- 
tions, except for the tt production mode for which only 
the NLO cross-section is available. The yields of expected 
SM Higgs signal as a function of the Higgs boson mass are 
reported in Table [T] The Higgs signal yield is dominated 
by gluon fusion while vector boson fusion contributes to 
about 10% of the total production cross-section. 
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Fig. 1. Invariant mass of the diphoton candidates for the 5063 
events passing the selection found in the 1.08 /fo"' of ATLAS 
data (black dots). The background model function (defined in 
Section |6f is represented by the red line. Finally, five times the 
SM Higgs signal expectation for mn = 120 GeV is represented 
by the dashed red line. 



4 Data sample composition 

The main background components to our study are the 
diphoton production, the photon-jet production for which 
the jet is misidentified as a photon, the dijet production 
with two jets misidentified as photons and finally the Drell- 
Yan continuum where both electrons are misidentidied as 
photons. 

To estimate the contributions of the events with at least 
one jet to the dataset, a method based on control regions 
of data was used. It relies on two discriminating variables: 
the isolation and the photon identification cut. The rejec- 
tion powers of these variables are mostly independent. The 
method itself is a generalization of the method used in Ref. 
||9 |. The Drell-Yan background is estimated from Z — > ee 
candidate events, by drawing the probability of misidenti- 
fying one of the two electrons as a photon. This probability 
is then multiplied to the Drell-Yan predicted yield to eval- 
uate the corresponding background rate. 
In the diphoton mass range 100 to 160 GeV, the number 
of true diphoton events is found to be 3650 + 100 + 290, 
where the first uncertainty is statistical and the second sys- 
tematic. The photon-jet and jet-jet contributions are esti- 
mated to be 1110 + 60 ± 270 and 220 + 20+ 130 events, 
respectively. The Drell-Yan background, which has a softer 
invariant mass spectrum than the other components, is es- 
timated to 86 ± 1 + 14 events. The different contributions to 
the diphoton mass range distribution are shown in Fig.|2] 

5 Event categorization 

Several topological and kinematical event categories are 
defined in order to benefit from the different signal over 
background ratios in the different regions of phase space. 
The present analysis relies on 5 coupled 77 region/conversion 
status categories: 

1 . Unconverted central : both photons unconverted and both 
located in the central region of the EM calorimeter, i.e. 
with It/I < 0.75; 

2. Unconverted rest: both photons unconverted but at least 
one is not located in the central region; 
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Fig. 2. Diphoton, photon-jet, dijet and Drell-Yan contributions to 
the diphoton candidate invariant mass distribution. The various 
components are stacked on top of each other. The error bars cor- 
respond to the uncertainties on each component separately. 

Table 1. Yields of Higgs signal at the MC generated mass points. 



m„ [GeV] 


110 


120 


130 


140 


150 


Signal yield 


17.0 


17.6 


15.8 


12.1 


7.7 



3. Converted central: at least one photon is converted but 
both are located in the central region; 

4. Converted transition: at least one photon is converted 
and at least one photon is located in the region where 
1.30 < M < 1.75 

5. Converted rest: at least one photon is converted and the 
event does not lie in the two previous categories. 



6 Bacl<ground and signai modeis 

The background distributions in each category are esti- 
mated by fitting the diphoton mass distribution by an ex- 
ponential function in the whole selected range. 
The MC Higgs signal diphoton mass distributions were 
fitted using a model consisting in a Crystal Ball, to de- 
scribe the core of the distribution, and a gaussian to de- 
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Fig. 3. Fit result of the Higgs signal diphoton invariant mass, for 
niM = 120 GeV, for all the categories together The core compo- 
nent of the mass resolution, (Tcb, is 1.7 GeV and the FWHM of 
the distribution is 4.0 GeV. 

Table 2. Systematic uncertainties on signal. 



Uncertainties on the predicted signal yield 

Photon reconstruction and identification (ID) 
Isolation cut efficiency 
Trigger efficiency 
Signal cross section 

Signal acceptance from modeling of the Higgs 
Luminosity 



+ 11% 
+3% 
±1% 

+20% 

-15 

±1% 

±3.7% 



Uncertainties on the mass resolution 




Calorimeter energy resolution 


±12% 


Energy calibration extrapolation from e to y 


±6% 


Effect of pileup on energy resolution 


±3% 


Photon angle measurement 


±1% 
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Fig. 4. The observed and expected 95% confidence level limits, 
normalized to the SM Higgs boson cross sections, as a function 
of the Higgs boson mass. 



The combined likelihood function is the product of the 
likelihood functions for the 5 categories. The systematic 
uncertainties are treated by using nuisance parameters fol- 
lowing gaussian PDFs. The 95% confidence level limit on 
of the inclusive production cross section of a SM-like Higgs 
boson relative to the SM cross section, is shown on Fig. 
|4] The expected CLs spans from 3.3 to 5.8 times the SM 
production cross section for masses between 110 and 150 
GeV. The observed limits are between 2.0 and 5.8 times 
the SM cross section. No significant excess is observed. 
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scribe the tails of the distribution. A fit to the Higgs bo- 
son mass shapes is performed on the distributions given 
by the Monte-Carlo samples. For the simulated point with 
niH = 120 GeV, the core components of the mass resolu- 
tion have a width which spans between 1.4 and 2.1 GeV, 
depending on the category. The fit result for the nin =120 
GeV sample is represented on Fig.|3]for the inclusive sam- 
ple after selection. 

7 Systematic Uncertainties 

As far as Higgs signal is concerned, two types of system- 
atic uncertainties were identified: uncertainties affecting 
the signal yield and uncertainties concerning the invariant 
mass resolution. They are summarized in Table |2] Photon 
identification and energy resolution systematics dominate. 
The background modeling uncertainties are estimated by 
checking how accurately the chosen model fits different 
predicted diphoton mass distributions and comparing dif- 
ferent functional forms for the background model. The re- 
sulting uncertainty is between +5 events at 110 GeV and 
±3 events at 150 GeV for a Higgs boson mass signal region 
about 4 GeV wide. 



8 Results 

The statistical approach used for this analysis is the mod- 
ified frequentist approach (CLs) documented in Ref. IfTOl . 
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